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ABSTRACT 

We present an algorithm for selecting an uniform sample of gravitationally lensed quasar candidates 
from low-redshift (0.6 < z < 2.2) quasars brighter than i = 19.1 that have been spectroscopically 
identified in the Sloan Digital Sky Survey (SDSS). Our algorithm uses morphological and color se- 
lections that are intended to identify small- and large-separation lenses, respectively. Our selection 
algorithm only relies on parameters that the SDSS standard image processing pipeline generates, al- 
lowing easy and fast selection of lens candidates. The algorithm has been tested against simulated 
SDSS images, which adopt distributions of field and quasar parameters taken from the real SDSS 
data as input. Furthermore, we take differential reddening into account. We find that our selection 
algorithm is almost complete down to separations of 1" and flux ratios of 10~ . The algorithm selects 
both double and quadruple lenses. At a separation of 2", doubles and quads are selected with similar 
completeness, and above (below) 2" the selection of quads is better (worse) than for doubles. Our 
morphological selection identifies a non-negligible fraction of single quasars: To remove these we fit 
images of candidates with a model of two point sources and reject those with unusually small image 
separations and/or large magnitude differences between the two point sources. We estimate the effi- 
ciency of our selection algorithm to be at least 8% at image separations smaller than 2", comparable 
to that of radio surveys. The efficiency declines as the image separation increases, because of larger 
contamination from stars. We also present the magnification factor of lensed images as a function of 
the image separation, which is needed for accurate computation of magnification bias. 
Subject headings: gravitational lensing — quasars: general 



1. INTRODUCTION 

Gravitational lensing is a unique tool for exploring cos- 
mology and the structures of astronomical objects. The 
well-understood physics of lensing makes it straightfor- 
ward to use as a cosmological probe. Moreover, it is the 
only method that probes distributions of dark matter di- 
rectly, since gravitational lensing is a purely gravitational 
phenomenon. 

While detailed investigations of single intriguing lens 
systems are useful for cosmological and astrophysical ap- 
plications, for some studies it is essential to do statis- 
tical analyses of a complete sample of lenses that are 
selected from a well-understood source population. For 
instance, strong lens probabilities are sensitive to the vol- 
ume element of the universe, and thus the amplitude of 
the cosmological constant, thus we can constrain cosmo- 
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logical models by comparing the observed lensed frac- 
tion in a catalog of distant sources to theoretical ex- 
pectations (lTurnerlll990HFukugita et alJll99fllKochaneld 



1199ft iChiba fc Yoshiil 11999ft . Strong lensing probabili 
ties at the cluster mass scale p robe the abundance and 
mass distribution of clusters ijNaravan fe White! 119881 
iMaoz et alJll997|) . The lens image separation distribu- 
tion from galaxy to cluster mass sc ales can be used to 
study baryon cooling in dark halos l)Kochanek fc White) 
12001ft and the connectio n between g alaxy luminosities 
and masses of host halos {Oeuri 200i|). 

A statistical sample is important not only for stud- 
ies of lensing rates, but also for interpretations of the 
results of individual lens modeling. It has been ar- 
gued that lensing objects may be more or less atyp i- 
cal in their shapes , orientations (Henna wi et alJ l2006b1. 
and environments l)Qguri et alJl2005b|) because of strong 
dependences of lensing probabilities on these quanti- 
ties. Therefore we must take proper account of lens- 
ing biases to make a fair comparison between the- 
ory and observed results, which can only be accom- 
plished for a well-defined statistical lens sample. A sim- 
ilar argument holds for constraints from stacked sam- 
ples of strong lenses, including the structure and evo- 
lution of early-type gala xies l)Treu fc Koopmansl 120041: 
iRusin fc Kochanekl 12005ft and the fraction of su bstruc- 
tures in lens galaxies ijDalal fc Kochanekl 12002'). It is 
quite hard to estimate the degree of the lensing biases if 
one adopts a heterogeneous sample of lenses that have 
been discovered in a variety of ways. 

The current largest complete sample of strong lenses 
was selected from a radio survey: The Cosmic-Lens All 
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Sky Survey (CLASS; IMvers et all l2003 IBrowne et all 
2003) has discovered 22 gravitational lenses from 
~16,000 radio sources. 13 of these lenses out of ^9,000 
flat-spectrum radio sources constitute a statistically well- 
defined sample. Recent progress on large-scale optical 
surveys has made it possible to construct comparable or 
even larger lens samples in the optical band. In par- 
ticular, the Sloan Digital Sky Survey fSDSS: lYork et all 
I2000D is expected to id entify ^100,000 confirmed quasars 
({Schneider et al|2005l) and ~500, 000 photometrically se- 
lected quasars ijRichards et al.l 12004) when completed. 
In addition to the large number of sources, the opti- 
cal lens survey has the advantage of a source popula- 
tion whose redshift and magnitude distributions are well 
understood, which is important when doing statistical 
analyses. On the other hand, optical lensed quasars are 
sometimes significantly affected by the light from fore- 
ground lens galaxies and dust reddening. Therefore, op- 
tical and radio lens surveys contain different systematics, 
and it is thus of great importance to compare the results 
from these two complementary lens samples. 

We are searching for optical gravitationally lensed 
quasars from the imaging and spectroscopic data of 
the SDSS to construct a statistical sample of lensed 
quasars. Indeed, this project, the SDSS Quasar Lens 
Search (SQLS), has already discovered 14 new lensed 
quasars, and has re covered several (~ 7) previously 
known lensed quasars flnada et alJ2003alblQl2005ll2006l: 
Johnston et a,l]l2003 IMorga.n et al.ll2003t IPindorTHh 
2004 1 1200ft IQgiiri et al.l I2004H I2005al IBurles et al. 
2006). Although these lenses do not constitute a statisti- 
cally well-defined sample, these discoveries give promise 
that we can create a large complete sample of lensed 
quasars. In order to construct a statistical sample, we 
need (i) to select lens candidates by a homogeneous 
method, and (ii) to understand the completeness of the 
method. In reality, we must also keep the efficiency rea- 
sonably high to assure the feasibility of follow-up ob- 
servations. We can take advantage of the deep SDSS 
imaging data in five broad bands fo r the efficient i denti- 
fication of lensed quasar candidates. Pin dor et al.l (2003, 
hereafter P03) presented a selection method for lensed 
quasars from the SDSS, and quantified its selection func- 
tion. Their selection algorithm compares the chi-square 
value of a two-component point-spread function (PSF) 
model of spectroscopically confirmed SDSS quasars with 
that of a one-component PSF model. From simulated 
SDSS images, they found that the algorithm recovers 
quasar pairs of similar fluxes down to separations of 
~ 0'/7. 

In this paper, we present another efficient algorithm to 
identify low-redshift (z < 2.2) lensed quasar candidates 
from the spectroscopic sample of SDSS quasars. Our 
selection algorithm uses both morphological and color 
selection. The former aims to select quasars that are ex- 
tended and are not well fitted by a PSF; it is designed 
to identify small-separation (~ 1") lens candidates. The 
latter algorithm simply compares colors of close pairs of 
objects with those of source quasars, and is focused on 
lenses whose image separations are large enough to be 
deblended in the SDSS imaging data reduction pipeline. 
Our selection algorithm relies only on the outputs of the 
SDSS imaging pipeline and therefore does not require in- 
dependent analysis of the imaging data themselves. This 



allows easy and fast identification of lensed quasar can- 
didates. We apply our algorithm to realistically simu- 
lated SDSS images to quantify the completeness in detail. 
For the simulation of lenses, we follow the methodology 
developed by P03. We adopt observed distributions of 
sky levels and seeing sizes, as well as the observed SDSS 
source quasar population in a real SDSS quasar catalog, 
to estimate the completeness of our algorithm in a real- 
istic manner. The algorithm will be applied to the SDSS 
data in Paper II (N. Inada et al., in preparation). 

This paper is organized as follows. In ^ we discuss 
the source quasar sample suitable for lens statistics. Our 
lens selection algorithms are presented in and are 
tested against simulations in Section [3] is devoted to 
study an additional selection that is intended to remove 
any point-source like candidates. The efficiency of our 
algorithms is briefly discussed in ^ and we summarize 
our results in f|7| 

2. SOURCE QUASAR SAMPLE 

The SDSS is a survey to image 10 4 deg 2 of 
the sky, and also to conduct spectroscopy of galax- 
ies, quasars, a nd stars that a r e sele ct ed from the 
imag i ng data ijEisenstein et al.1 I200H IStrauss et alJ 
120021 IRichards et alJ l2002j) . The tiling algorithm 
(Blanton et al. 2003) assign s targets to fiber s. 10 A ded- 
icated 2.5-meter telescope ijGunn et al.ll2006T) at Apache 
Point Observ atory (APO) is e quipped with a multi- 
CCD camera (|Gunn et al.lll998j) with five optical broad 
bands centered at 3551, 4686, 6166, 7480 and 8932 A 
l|Fukugita et al.lll99fil:fStoughton et all2002j) . The imag- 
ing data are auto matically reduced b y the photometric 
pipeline PHOTO l|Lupton et alJl200lD . The astrometric 
positions are acc urate to about O'/ l for sources brighter 
than r = 20.5 (|Pier et all 12003(1 . and the photomet- 
ric cali bration errors are typically less than 0.03 mag- 
nitude l|Hogg et al.ll200H: iSmith et a.l.ll2002t llvezic et alJ 
l20f)l. All da,ta are being released on a regular sched- 
ule (Stoughton et alJl2002t lAbazaiian et al.ll2003l 12004 
l2005HAdelman-McCarthv et alJl200^T 

We construct a statistical sample of lensed 
quas ars from the spectroscopic SPSS quasar cata- 
log l|Schneider eTaTl 120021 I20M 120051 . The quasar 
sample is quite heterogeneous in the sense that the 
quasars are selected with several different methods and 
the completeness is strongly dependent on the properties 
of quasars (e.g., magnitu des, redshifts, morphologies; 
see IRichards et al"l l2006bj) . We make a subsample of 
source quasars, which are suitable for lens statistics, by 
posing the following restrictions. 

1. We consider only quasars in the redshift range 
0.6 < z < 2.2. Lower-redshift quasars are often 
extended in the SDSS images, which makes the 
morphological selection fi J3.1|l quite difficult; in any 
case, these quasars have low lensing probabilities. 
At low redshift most quasars are targeted by UV 
excess which can select both point and extended 
quasars efficiently, while at higher redshifts outliers 
from the griz stellar locus are targeted as quasar 

10 Because of the physical size of fibers, spectra of no two targets 
that are closer than 55" are taken on a single plate. This "fiber 
collision" limitation makes it unlikely to obtain spectra of more 
than one component in a gravitational lens. 
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Fig. 1. — The completeness of the quasar spectroscopic target 
selection for poi nt and extended quasars , computed from the sim- 
ulation done by Richards ct al. (2006b). We show completeness 
averaged over 15 < i < 19.1. The ratio of the two completenesses 
is shown in the lower panel. 

candidates, for which we had to restrict our atten- 
tion to point sources because of large contamina- 
tion from galaxies. This is clearly shown in Figure 
^ which plots the completene s s obtai ned from the 
simulation by iRichards et alJ l)2006b|) . The com- 
pletenesses are different at z > 2.2 between point 
and extended quasars, which introduces a non- 
negligible anti-lens bias. In particular, extended 
quasars at z > 3 are not selected by the quasar 
target selection algorithm, which significantly re- 
duces the number of ~ 1" separation high-redshift 
lenses in the SDSS spectroscopic quasar catalog 
(see also Appendix El . In addition, at z > 2.2 
it becomes more difficult to use colors to distin- 
guish quasars from stars, thus the efficiency of lens 
candidate selections is significantly degraded, par- 
ticularly at the large-separation region where con- 
tamination from stars are quite large. 

2. We ado pt quasars with Ga lactic extinction cor- 
rected llSchlegel et alJ 1199811 i-band magnitudes 
15.0 < i cor < 19.1. This is because at i COI < 19.1 
the low-redshift [z < 2.2) quasar target selection 
is almost complete ^Richards et all 120021 H)06b), 
while fainter quasars are select ed mainly with 
vario us "serendipity" algorithms ijStoughton et al.1 
2002) which are very incomplete and may introduce 
biases in estimating lensing rates. Moreover, the 
faint quasars have several disadvantages, including 
large photometric errors that reduce the efficiency 
of lens selection algorithms, and contamination by 
lens galaxies that beco me more importa nt for less 
luminous quasars (e.g., lKocha^el3[T99^ . Fainter 



lenses are identified as lens candidates with lower 
efficiency because of their lower signal-to-noise ra- 
tios (see • 

3. Close lens pairs are difficult to distinguish from sin- 
gle PSFs if the seeing is poor. Therefore we reject 
quasars in fields if the i-band image has the SDSS 
imaging parameter PSF_WIDTH (effective PSF width 
based on 2-Gaussian fit) > l'/8. This criterion re- 
jects roughly 2% of quasars. 

As a specific example, we adopt the 46,420 quasars 
in the DR3 quasar sample ({Schneider et al.ll2005|) as our 
parent sample. By restricting the range of i-band mag- 
nitudes and redshifts, the number of quasars decreases 
to 23,316 (~ 50%). The final number of quasars, after 
rejecting poor seeing fields, is 22,682 (~ 49%). 

3. SELECTION ALGORITHM 

In this section, we present our two selection algorithms, 
using morphological and color selection criteria, respec- 
tively. The algorithms are intended to identify low- 
redshift (z < 2.2) lensed quasars, though their extension 
to higher-redshift lenses is rather straightforward (be- 
sides the inefficiency and imcompleteness of high-redshift 
lenses). Lens candidates are identified using imaging 
data; we do not use any spectroscopic information, al- 
though in some cases spectra off er a powerful way to 
locate gravitational lens ing (e.s.. Uohnston et alJ 120031: 
iBolton et all I2005L 120061) . We determine the selection 
criteria described below mostly in an empirical manner; 
we choose criteria to keep the completeness high, for both 
known lenses and simulations done in 21 and at the same 
time have reasonable efficiency when applied to the SDSS 
data. 

3.1. Morphological Selection 

When the image separation between multiple images 
is small (6 < 2'.'5), PHOTO cannot deblend the system 
into two components, and the objects are classified as sin- 
gle extended objects. Thus small-separation lenses can 
be identified as lens candidates by searching for objects 
that are poorly fitted with the local PSF. To do so, we 
use the following two SDSS imaging parameters. One is 
obj c_type which describes the classification of the object 
as stars or galaxies, determined from the differences of 
magnitudes obtained by fitting PSF and galaxy profile s 
to the images in each band fe.g.. lStoughton et alJ l2002): 
objc_type = 6 indicates the object is a point source, 
while objc_type = 3 means the object is extended. The 
other parameter, star_L (available in each band), is the 
probability that an object would have at least the mea- 
sured value of x 2 of a fit of the image to the PSF, if it 
really is well represented by a PSF. Put another way, the 
higher the value of star_L, the more likely it is fitted by 
a PSF. 

The specific selection criteria of this morphology selec- 
tion are as follows. First, even if a quasar is classified 
as a point source, we select it as a lens candidate if the 
PSF likelihood is small. Specifically we select objects 
that satisfy all of the following four conditions: 

(Ml): objc.type = 6, 

star_L(it) < 0.03, 
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star_L(g) < 0.04, 

star_L(r) < 0.07 OR star_L(i) < 0.07. (1) 

We do not use the z-band star_L parameter because of 
the relatively low signal-to- noise ratio (S/N) of z-band 
images. Although the S/N of u-band images are also 
low, we use the it-band parameter because quasars are 
UV-excess sources and thus the M-band is quite effective 
in discriminating quasar-quasar pairs from quasar-star 
pairs. However, we find that the condition on the it-band 
parameter is sometimes too strong, resulting in missing 
some lens candidates. As a separate cut, we adopt a re- 
laxed range of u-band star likelihood, but slightly tighter 
conditions on other parameters: 

(M2): objcjtype = 6, 

star_L(it) < 0.06, 
star_L(.g) < 0.04, 

star_L(r) < 0.04 OR star_L(i) < 0.04. (2) 

When a quasar image is classified as extended, we use 
a similar but relaxed, criterion to identify a lens candi- 
date; an object is selected when all the following three 
conditions are met: 

(M3): objc_type = 3, 

star_L(w) < 0.45 OR star_L(.g) < 0.35, 
star_L(Y) < 0.60 OR star_L(i) < 0.60. (3) 

We select quasars that satisfy one or more of the crite- 
ria (M1)-(M3) as lens candidates. In practice, approx- 
imately half of small-separation gravitational lenses are 
classified as extended objects and selected by the con- 
dition (M3). Lenses that are classified as point-sources 
tend to satisfy both conditions (Ml) and (M2); the differ- 
ence between (Ml) and (M2) becomes important only for 
limiting cases in which lenses are close to the magnitude 
limit (i ~ 19.1) and/or have very small image separa- 
tions (6 < 1"). These conditions are designed using both 
real SDSS data and simulations presented in 21 

3.2. Color Selection 

When the image separation is large enough for PHOTO 
to deblend two components, we can select lens candidates 
by comparing the colors between spectroscopically con- 
firmed quasars and candidate companions. At 6 > 7" 
the candidates' companions must be point sources, but 
for smaller separations we also allow extended sources to 
be candidate companions, because for such small separa- 
tions lensing galaxies may be superposed on the fainter 
components (e.g., Q 0957+561). The selection method 
presented here is s ome what similar to those given by 
lOeuri et all l)2004aj) and lHennawi et alJ l)2006aj) . but we 
follow P03 to extend the methods to select differentially 
reddened lenses. 

We conduct the color selection on the basis of the color 
difference: 

D(i _y) — (z j)quasar (^ j)companion 

OR (i j)companion (^ J/quasarj (4) 

where (i, j)=(u, g), (g, r), (r, i), (i, z), and quasar 
and companion indicate the SDSS spectroscopic quasar 
in the source quasar sample and candidate companion, 
respectively. We use PSF magnitudes and their errors 



throughout the paper. Although the definition of D(i— j) 
contains a sign ambiguity, we regard candidate pairs as 
lens candidates if they pass the criteria below using either 
one of the above definitions: This becomes important 
when we consider differential reddening, because at first 
sight we do not know which object is reddened more. 
When the color difference comes from reddening, D(i—j) 
for all set of (i, j) must have the same sign. On the other 
hand, the color difference of two unrelated objects could 
have both positive and negative signs for different (i, 
j). This is the reason we do not define D(i — j) as the 
absolute value of the color difference. 

Since gravitational lensing does not change the color 
of objects, the color difference should be small if two 
components are lensed images of a single quasar. This 
criterion can be written as follows: 

(CI): \D(i — j)\ < AD(i — j), (5) 

which must be met for all four colors. The limit AD(i — 
j) describes an acceptable error of the color difference, 
which we adopt the following fixed values 



AD(u 


-9) 


= 3 x 0.20, 


AD(g 


-r) 


= 3 x 0.09, 


AD(r 


-i) 


= 3 x 0.09, 


AD(i 


-*) 


= 3 x 0.30. 



(6) 



These values correspond to typical 3-cr errors of the color 
differences of faint (~ 19 — 20 mag) objects in the SDSS 
images. We find that color differences contain very large 
errors when image separations are small because of the 
SDSS deblending algorithm. To account for this, we re- 
place AD(i - j) with 3 x (TD{%-j) f° r 2 " < 6 < 3'/5, and 
with 6 x cr£)(j_j) for 8 < 2", where arj(i-j) is defined by 

°~D(i—j) = \J (of + Cr|)quasar + (°f + fj) companion i (7) 

where o~i are PSF magnitude errors estimated by 
PHOTO. Since it is quite rare that objects lie at so close 
to the quasar, the large errors have little effect on the 
efficiency of our lens selection algorithm. 

While gravitational lensing is independent of wave- 
length, significant color differences between lensed com- 
ponents may be caused by reddening. For instance, 
iFalco et al.l l)1999[) measured the distribution of differ- 
ential reddening AE(B — V) from a number of gravi- 
tationally lensed quasar systems and found that is well 
described by a Gaussian with a zero mean and standard 
deviation of 0.1 mag. We estimate how the differential 
reddening affects the color difference D(i — j) as follows 
(see also P03). First, we use the co mposite quasar spec- 
trum of IVanden Berk et all 1)2001^ as an input quasar 
spectral energy distri bution. Next we ado pt the Galac- 
tic extinction law 11 of lCardelli et al.l l)1989T) and compute 
the color difference as a function of quasar and lens red- 
shifts and AE(B — V), assuming Ry = 3.1. We show 
estimated reddening vectors for several sets of quasar and 
lens redshifts in Figure [3 As expected, the effect of the 
differential reddening is more significant for bluer bands. 

11 We note that our result is not strongly dependent on the as- 
sumed reddening law. For instance, we have repeated the a nalyses 
adopting SMC-like reddening, which [Hopkins ct al. 12004) claims 
fits the observed reddeni ng of quasars, and confirmed that the dif- 
ference of completeness f t|4.3i is indeed small. 
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Fig. 2. — Reddening vectors in color difference-color difference space are shown by solid lines plus filled triangles. We compute the 
vectors for the quasar redshift z s = 0.7, 1.2, 1.7, 2.2, and the lens redshift z s = 0.2, 0.4, 0.6, 0.8. The differential extinction AE(B - V) 
is changed from to 0.3 (triangles are plotted every 0.1 increment). The color selections (CI) and (C2) are shown by dashed and dotted 
boxes, respectively. 



From this estimation, we prepare another color selection 
algorithm in color-color space that is designed to select 
lens systems with large differential reddening: 
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(9) 



xi + Ax 

x < 0.2 + Aa; IF x = D(u - g), (8) 

where the condition must be satisfied for all set of (x, 
y) = {D{u-g),D{g-r)), (D(g-r), D(r-i)), (D(r-i), 
D(i — z)), and 

{x 1 ,y 1 ,x 2 ,y2) 
( (0.3, 0.4, 0.4, 0.4) (x, y) = (D(u ~ g), D(g - r)), 
= I (0.2, 0.2, 0.4, 0.1) (x, y) = (D(g - r),D(r - i)), 
{ (0.1, 0.2, 0.4, 0.1) (x, y) = (D(r - i), D(i - z)). 

The color difference errors, Aa; and Ay, were defined 
in equation © and the text that follows it. Figure |21 
shows the color selection regions defined by (CI) and 
(C2). In particular, Figure indicates that the condition 
(C2) can identify lens systems with very large differential 
reddening, AE(B - V) ~ 0.3. The cut at D(u - g) = 
0.2 + AD(u — g) is added to reduce the contamination of 
quasar-star pairs. We select quasars that satisfy one or 
both the criteria (CI) or (C2) as lens candidates. 

4. TESTING SELECTION ALGORITHM 
4.1. Simulation Method 



To quantify the completeness of our selection algo- 
rithm, we perform a series of simulations that mimic 
SDSS observations. We follow the methodology devel- 
oped by P03, but in this paper we implement multicolor 
(ugriz) simulations to test our selection based both on 
the morphology and color match. 

First, we prepare object fields in which simulated 
lensed images are placed. We assign the characteris- 
tics of each field (sky level, seeing, and magnitude zero- 
point) from a randomly selected field in the real SDSS 
data. Figure |31 shows the sky counts and seeing sizes 
(PSF.WIDTH) of 1000 random fields. Following the dis- 
cussion in we reject fields with PSF_WIDTH > 1'.'8 in 
i-band in our simulations. 

For each field, we also choose a quasar randomly from 
the subset of the DR3 quasar catalog with 0.6 < z < 2.2 
and 15.0 < i cor < 19.1, and adopt its redshift and magni- 
tudes in the simulation. However, the redshift and mag- 
nitude distributions of the lensed population are quite 
different from those of the original unlensed population, 
because the lensing probability is a strong funct ion of 
redshift and mag nitude (e.g.. iTurner et ali n"9841. We 
account for this difference by assigning to each object 
a weight proportional to its lensing probability, com- 
puted from a singular isothermal sphere model for the 
mass distribution of lensing galaxies and the velocity 
function of early-type galaxies derived from the SDSS 
l|Sheth et al.ll200l iMitchell et all 12005]) . The velocity 
function in comoving units is assumed to be constant 
with redshift. The magnification bias is computed by 
assuming the luminosity function of quasars obtained 
from th e combination of the T wo-Degree Field and SDSS 
results (|Blchards et alJ 120051 "2SLAQ+Croom et al." 
model). Although the lensing probability is sensitive to 
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Fig. 3. — The width of PDFs versus sky levels for 1000 randomly selected fields in the SDSS for each band. Open circles indicate those 
fields that are rejected because of poor seeing. 
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Fig. 4. — Rcdshifts versus i-band magnitudes of randomly se- 
lected quasars at 0.6 < z < 2.2 and 15.0 < i CO r < 19-1 are plot- 
ted. Open circles are selected from the quasar sample without any 
weight, and filled triangles used a lensing probability weight. 

cosmological parameters, their main effect is to change 
the overall amplitude of the lensing probability and to 
have little effect on the dependence of lensing probabil- 
ity on the source redshift and magnitude, on which we 
are focusing. Figure^lplots rcdshifts versus «-band mag- 
nitudes for randomly selected quasars with and without 
lensing probability weighting. Because of magnification 
bias, gravitational lensing preferentially chooses brighter 
and higher-redshift quasars. 

We convert magnitud e to the number of counts with 
the following equations l|Lupton et al.lll999(h 



2.5 
"ln(10) 

F _ 



asinh 



{F/Fo 



V 2b 
counts 



ln(&) 



exptime|sec| 



1 Q0.4(Z+fcX) 



(10) 
(11) 



where b is the softening parameter, Z is the magni- 
tude zero-point, k is the extinction coefficient, X is air- 
mass, and the exposure time is 53.91 sec. Gravitation- 
ally lensed images are created in the simulated image by 
placing two point-source objects that are separated by 
an angle 9 and have a flux ratio fi. For the point-spread 
function (PSF), we use an analytic PSF composed of 
two Moffat functio ns with j3\ = 7 and = 2 and widths 
given for each field <|Racindll996t P03). 12 The PSF mag- 
nitude of each system is normalized to that of the input 
quasar PSF magnitude in each band. With this proce- 
dure, we can simulate lensed quasars with realistic colors 
and field conditions. 

We also include differential reddening in our simula- 
tion. The value of differential extinction AE(B — V) 
is randomly distribut ed as a Gaussian with a standard 
deviation of 0.1 mag l)Falco et al.ll 1999]) . Again, the red- 
dening vecto rs are computed from the composite quasar 
spectrum of IVanden Berk et alJ (I2001T) and the extinc- 
tion curve of iCardelli et alJ (1989) assuming Ry = 3.1. 
The lens redshifts are also assigned randomly in the range 
[0.2, z max ] with z max — min(0.8, z s — 0.2), where z s is the 
redshift of the source quasar. The distribution has the 
average redshift z = 0.5, which is consistent with the 
value z ~ 0.51, the mean of 30 existing gravitationally 
lensed quasar systems with known lens redshifts. We 
note the flux ratio fi depends on the band when redden- 
ing is included; hereafter we refer to fi as the flux ratio 
of the z-band image. 

We next pass these simulated fields through the 
PHOTO pipeline. Here we briefly summarize the pro- 
cedure; for more details, see P03. We begin by prepar- 
ing a PSF frame in which bright PSF stars (constructed 

12 The width of the PSF measured by PHOTO (PSF.WIDTH) is 
systematically larger than than the input FWHM of the Moffat 
function by a factor 1.15 (P03). We take this difference into account 
throughout the paper. 
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Fig. 5. — The distributions of star likelihood star_L for lensed and single quasars, estimated from our simulations. For lensed quasars, 
we consider only those with objc.type = 6, image separations 9 = 1" — 1.4", and flux ratios /; > 10 -05 . Vertical dashed lines indicate 
the limits used in our morphological selection (Ml) and (M2) (see H3.ll . 



from two Moffat functions) are displayed. This PSF field 
is passed through the fake stamp pipeline and postage 
stamp pipeline to generate a psField file in which the re- 
sult of the PSF fit is stored. The object field is passed 
through the fake stamp pipeline and frames pipeline; 
this produces tables of the measured object parameters 
(fpObjc) as well as the cropped atlas images (fpAtlas). 
We do not need to pass the result through the photomet- 
ric calibration pipeline, since we can convert counts to 
magnitudes again by using equations 1)10)1 and 1)11(1 . The 
magnitude error is estimated from that of the counts as 



Am = 



2.5 A(counts) 1 ipQ-Hz+kX) 
ln(10) exptime[sec] 2b 



1 + [(F/F )/2bY 



(12) 

The image parameters used to select the lens candidates 
are extracted from the fpObjc file. When PHOTO de- 
blends two images, we regard the brighter image as tar- 
geted by the quasar spectroscopic target selection, and 
use the image parameters of the brighter image. 

We repeat the simulation for different sets of field pa- 
rameters, redshifts and magnitudes of quasars, and red- 
dening, and compute the completeness of our selection 
algorithm as a function of the image separation 9 and 
flux ratio Each simulation run consists of simulations 
of 25 x 21 quasars; the separation is changed from 0" to 
4'.'8 with a step of 0'.'2 and the flux ratio from lO^ 1 to 1 
with a step of a factor of 10 05 . Since the maximum im- 
age separation 9 = 4'.'8 is large enough compare with the 
seeing sizes of the SDSS images, we think the complete- 
nesses of larger-separation lenses can be approximated 
by those of 9 = 4'.'8. We note that only a few known 



lensed quasars have the flux ratios fi < 0.1. We sim- 
ulate 1000 runs to make an accurate estimation of the 
completeness. 

We note that our simulations differ from those of P03 
in several ways. In P03, the selection algorithm was 
tested for a few choices of seeing sizes and magnitudes 
of quasars, while in this paper we distributed them in a 
realistic way. We also simulated images in all five bands, 
while P03 did not consider differences between different 
bands. 

4.2. Star Likelihood 

Before showing the completeness, we check the distri- 
bution of star_L to see the validity of using this parame- 
ter to identify lens candidates. All single quasars should 
be classified as point sources (objc_type = 6), thus we 
can choose quite loose conditions for selecting lens can- 
didates, without a significant decrease of efficiency, when 
they are classified as extended objects (objc_type = 3; 
see (M3) in i )3.1fl . Therefore, in this section we concen- 
trate our attention on lenses with ob j c_type = 6, for 
which we need to adopt a good discriminator between 
lenses and single quasars. 

In Figure El we show the distributions of star_L in it, 
<7, r, and i for both lensed and single quasars from our 
simulations. For lensed quasars, we only consider the im- 
age separation range of 9 = 1" — 1.4" and the flux ratio 
fi > 10~ ' 5 . It is clear that the distribution of star_L for 
lensed quasars is markedly different from that for single 
quasars. We find that lensed quasars have quite small 
star likelihoods star_L ~ even when lenses are classi- 
fied as point sources. On the other hand, the distribution 
of star likelihoods for single quasars peaks at star_L ~ 1 
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0[arcsec] 

Fig. 6. — The completeness of our lens selection algorithm in 
the image separation fl-flu x ratio fi plane, computed from the sim- 
ulation described in i|4.1l From thick to thin lines, we plot com- 
pleteness contours of 0.95, 0.9, 0.75, 0.5, 0.25, 0.1, and 0.05. The 
dotted lines indicates the limit above which the selection is almost 
complete. 

and rapidly decreases as star _L decreases. This Figure 
assures that the conditions (Ml) and (M2) presented in 
H3.ll serve as an excellent method to distinguish gravita- 
tional lenses from normal quasars. 

4.3. Completeness 

We show the completeness in the 9-fi plane in Figure 
|SJ As expected, lenses with larger image separations and 
flux ratios closer to unity are more likely to be selected. 
In particular, our selection algorithm is > 95% complete 
at 6 > 1" and fa > lO" 5 (indicated by dotted lines). 
Therefore, we propose to construct a statistical sample 
within this region of parameter space. 

Figure shows the completenesses for three differ- 
ent ranges of redshifts/magnitudes. We find that the 
completeness is almost independent of the quasar red- 
shift. On the other hand, the completeness shows a weak 
but significant trend with magnitude, in the sense that 
brighter quasars are more complete. In particular, the 
completeness at small fi is larger for brighter quasars. 
This is reasonable since companion pairs, which are of- 
ten hidden when fi is small, are more easily detected 
for brighter quasars. However, even for faintest quasars 
(i > 18.7) the completeness at 9 > 1" and f { > I(T - 5 is 
very high. Therefore we conclude that our selection does 
not affect the distributions of redshifts and magnitudes 
of lensed quasars very much. 

Figure El does not make it clear if lenses are mainly 
chosen by the morphological or color selection at each 
point in the 9-fi plane. Figure [3] shows the completeness 
from the two selection algorithms separately. We also 
show the overlap of these two selections. At 9 < 2" 
lenses are mostly selected morphologically, and at 9 > 
2" the color selection dominates. The overlap of these 
two selections is quite small. The completeness from the 
two selection algorithms as well as total completeness are 
listed in Table □ 

4.4. Magnification Factor 
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TABLE 1 

Completeness for double lenses 



log(/i) 


e 

[arcsec] 




Completeness 




total 


morph 


color 


both 


— 1.00 


0.0 


0.036 


U.UoO 


U.UUU 


0.000 


-1.00 


0.2 


0.031 


0.031 


0.000 


0.000 


-1.00 


0.4 


0.032 


0.032 


0.000 


0.000 


-1.00 


0.6 


0.041 


0.041 


0.000 


0.000 


-1.00 


0.8 


0.100 


0.100 


0.000 


0.000 


-1.00 


1.0 


0.213 


0.213 


0.000 


0.000 


-1.00 


1.2 


0.304 


0.304 


0.000 


0.000 


-1.00 


1.4 


0.378 


0.378 


0.000 


0.000 


-1.00 


1.6 


0.415 


0.415 


0.001 


0.001 


-1.00 


1.8 


0.461 


0.459 


0.005 


0.003 


-1.00 


2.0 


0.472 


0.463 


0.011 


0.002 


Note. 


- Full table to 


appear in 


the on- 


line edi- 



tion. 

The magnification bias is one of the most important 
elements in predicting l ensing probabilities. The SDSS 
quasar target selection (Richar ds et alJ 120021) uses PSF 
magnitudes for the magnitude limit, therefore we exam- 
ine those of simulated lensed quasars with various im- 
age separations. Basically, when the image separation 
is quite large and two images are deblended successfully 
by PHOTO, the PSF magnitude of the targeted quasar 
should be roughly equal to that of the brighter image. 
However, when the image separation is too small to be 
distinguished from a point source, the PSF magnitude is 
that of the sum of the two images. Bearing this in mind, 
we define the following parameter: 



where /ihri and /Xf a ; are the magnification factors for the 
brighter and fainter images, respectively. The magnifi- 
cation factor fj, indicates the one measured in the SDSS 
(i.e., the PSF magnitude of the lens system divided by 
the original magnitude of the source quasar) . The defi- 
nition indicates that \i = [itot = /ibri + Mfai and /j = //bri 
correspond to p = 1 and p — 0, respectively. Therefore, 
the parameter should approach zero at large separations, 
and p — > 1 for 9 — > 0. 

We derive p as a function of the image separation from 
our simulation runs. We only consider simulated lens 
systems with flux ratio fi > 10~ 05 , because the derived 
p contains large errors when flux ratios are very large 
jt/bri 3> /Ufai > an d also because the completeness is low for 
large flux ratio systems. Moreover, /x to t ~ Mbri at /j <C 1 
implies that the magnification factors do not change as 
a function of the image separation in any case. Figure 
plots the result. As expected, p changes from p ~ 1 to 
~ with increasing image separations. We find that the 
curve is fitted roughly by the following form: 

p = - [1 + tanh (1.76 - 1.780)] , (14) 

where the image separation 9 is in units of arcsec. 

We note that for a given image separation the magni- 
fication bias can be written as 

B = pB tot + (l-p)B hri , (15) 

where B to t and -Bbri are the magnification bias from the 
total magnification and magnification of the brighter im- 
age, respectively. 
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Fig. 7. — The completeness of our lens selection algorithm for different redshift and magnitude bins. The definition of lines is same as 
Figure lUl Upper panels: The completenesses for restricted redshift ranges of z < 1.5 (left), 1.5 < z < 1.8 (center), and 1.8 < z (right). 
Lower panels: The completenesses for restricted i-band magnitude ranges of i < 18.2 (left), 18.2 < i < 18.7 (center), and 18.7 < i (right). 
In both cases, the sizes of bins are determined such that the (lensing probability weighted) number of quasars in each bin becomes roughly 
equal. 
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Fig. 8. — Left: The completeness of the morphological selection only. Center: The completeness of the color selection only. Right: The 
fraction of lenses selected from both morphological and color selections. The definition of lines is same as Figure 



4.5. Quadruple Lens Systems 

Thus far we have concentrated our attention on two- 
component lenses, but in reality a non-negligible fraction 
of gravitationally lensed quasars have four components. 
In this subsection we consider our selection algorithm's 
effectiveness for such quadruple lenses. 

The difficulty in studying the selection of quadruple 
lenses lies in their complexity; while double lenses are 
parameterized by flux ratio and image separation, we 
require 8 parameters to characterize a quadruple lens, 
making a full exploration almost impossible. Instead, in 
this paper we adopt the image configurations and flux 



ratios from 20 known quadruple lenses taken from the 
CASTLES webpage (see Table |2J). We apply a similar- 
ity transformation to the quadruple lenses in order to 
derive the completeness as a function of the image sep- 
aration. We define the image separation to be the max- 
imum distance between multiple images. The differen- 
tial reddening is in cluded, again using the distribution of 
iFalco et alJ ((1299). As in the case of double lenses, we 
simulate 1000 runs to make an accurate estimate of the 
completeness. 

Figure plots the completeness of each of the 20 
quadruple lens configurations as a function of the im- 
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TABLE 2 

QUADRUPLY LENSED QUASARS USED IN THE SIMULATION 



Name 



Type 3 



/qu 



quad 



Reference" 1 



PG 1115+080 
SDSS J0924+0129 
SDSS J1004+4112 
WFI2026-4536 
WFI2033-4723 
HE0230-2130 
MG0414+0534 
B0712+472 
B1555+375 
B1608+656 
H1413+117 
HE0435-1223 
HST12531-2914 
HST14113+5211 
HST14176+5226 

Q2237+030 
RX J0911+0551 
RX J1131-1231 
B0128+437 
B1422+231 



fold 2'.'43 0.17 CASTLES C 

fold l'/81 0.47 Keeton et al. (20061 

fold 14'.' 6 0.23 Oguri et al"72004al 

fold l'/44 0.31 Morgan et al. (2004') 

fold 2'.'53 0.56 Morgan et al. (2004*) 

fold 2'.'19 0.59 CASTLES 6 

fold 2'.' 13 0.47 CASTLES C 

fold 1'.'29 0.20 CASTLES C 

fold 0'.'42 0.49 Marlow et al. (19991 

fold 2'.'10 0.49 Fassnacht et al. (20021 

cross 1'.'35 0.71 CASTLES 6 

cross 2/56 0.62 Wisotzki et al. (20021 

cross 1'.'38 0.92 CASTLES 6 

cross 2'.'26 0.82 CASTLES 6 

cross 3'.'26 0.79 CASTLES 6 

cross 1'.'83 0.31 CASTLES 6 

cusp 3'.'25 0.31 CASTLES 6 

cusp 3'.'24 0.09 Sluse et al. (20031 

cusp 0'.'55 0.53 Phillips et al. (20001 

cusp 1'.'29 0.03 CASTLES 6 



a Based on visual inspection. 

b The maximum separation between multiple images. 

6 The flux ratio between the brightest image and the image 
farthest from the brightest image. 

d The reference from which we adopt the image positions 
and fluxes for the simulation. For fluxes, we use those of the 
j-band (F814W for the Hubble Space Telescope data) image 
when the optical data are available. If not, we adopt those 
of the radio band image. 

6 C. S. Kochanek, E. E. Falco, C. Impey, J. 
Lehar, B. McLeod, H.-W. Rix. CASTLES survey, 
|http://cfa" www.harvard.edu/castles/| 
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Fig. 9. — Th e m agnification factor of lensed quasars. The pa- 
rameter p, (eq. 1131 1 is plotted as a function of image separation 
0. For each image separation the mean and its standard d evia tion 
are shown. The dashed line indicates the analytic fit (eq. 1141 1. 

age separation, as well as the average of them. We find 
that our algorithm can select quadruple lenses reason- 
ably well: Even in the limiting case, 9 = 1", the com- 
pleteness is 82% on average, and at 9 > 2" it is al- 
most unity. We note that in both cases the complete- 
nesses vary most rapidly at 9 ~ 1" where the number 
of lensed quasars peaks. There is a systematic differ- 
ence of the completeness between double and quadruple 
lenses: At 9 < 2" the completeness for double lenses is 
larger than that for quadruple lenses, while at 9 > 2" 
the quadruple lenses are more complete. The differ- 




quad 

— quad, average 

- double (fjMO" - 5 ) 

L 



J I I I I I I I I I I I I I I I I I I I I I L 



12 3 4 

0[arcsec] 

Fig. 10. — The completeness for quadruple lenses as a function 
of the image separation. The dotted line shows the completeness 
of each quadruple lens image configuration and flux ratios taken 
from real observational data (see Table |5J . The thick solid line 
indicates the completeness averaged over the 20 dotted lines. The 
completeness of double lenses, averaged over the flux ratios in the 
range 10~ 0,5 < /, < 1, is also shown for reference. 

ences are small, but they might become important when 
discuss ing the relative number of quadruple and double 
lenses lIKeeton et alJll997tlRusin fc Teemarkl200lt IChad 
l27MlQgu7i fc Keetonll2004D . 

Figure ^2 shows the completeness for morphological 
and color selections separately. Unlike double lenses, 
quadruple lenses are effectively chosen by the morpho- 
logical selection even at large image separations, 9 > 3". 
In addition, it is quite common that lenses are selected by 
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Fig. 12. — The completeness for quadruple lenses of different 
image configuration types: fold (solid) , cross (dotted), and cusp 
( dashed) . 

both criteria This is because quadruple lenses are more 
complex than double lenses. In particular, the complex- 
ity and larger number of images for quadruple lenses 
make the completeness at 8 > 2" quite large. Again, 
the completeness from the two selection algorithms as 
well as total completeness are listed in Table 

We also study the completeness for different image con- 
figuration types: fold (two of the images lie close to- 
gether), cross (almost symmetric), and cusp (three of 
the images lie close together). To do so, we divide the 20 
quadruple lenses into these three categories by visual in- 
spection, and compute the average completeness for each 
category. The categorization is summarized in Table |2 
and the results are shown in Figure 1121 Our algorithm 
selects cross lenses very well, but is less effective for cusp 
lenses; this is because the fourth image of a cusp lens, 
which lies far from the critical curve, is usually much 
fainter than the other images near the critical curve and 
therefore the effective "flux ratio" of cusp lenses becomes 
small. To pursue this further, we consider the flux ratio 
between the brightest image and the image farthest from 
the brightest image (denoted by / qU ad), and see its corre- 
lation with the completeness. Figure IT51 plots / qua d and 
the completeness at 9 = 1" for all 20 quadruple lenses 
listed in Table As expected, we see a clear correla- 
tion: All quadruple lenses with completeness < 0.8 have 
/quad < 0.4. Therefore, despite the complexity and large 
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TABLE 3 

Completeness for quadruple lenses 



Name 



Completeness 







[arcsec] 


total 


morph 


color 


both 


PG1115- 


h080 


0.0 


0.035 


0.035 


0.000 


0.000 


PG1115- 


h-080 


0.2 


0.030 


0.030 


0.000 


0.000 


PG1115- 


1-080 


0.4 


0.029 


0.029 


0.000 


0.000 


PG1115- 


t-080 


0.6 


0.059 


0.059 


0.000 


0.000 


PG1115- 


t-080 


0.8 


0.188 


0.188 


0.000 


0.000 


PG1115- 


h080 


1.0 


0.490 


0.490 


0.000 


0.000 


PG1115- 


h080 


1.2 


0.796 


0.796 


0.000 


0.000 


PG1115- 


h-080 


1.4 


0.937 


0.937 


0.000 


0.000 


PG1115- 


t-080 


1.6 


0.978 


0.978 


0.001 


0.001 


PG1115- 


t-080 


1.8 


0.987 


0.987 


0.007 


0.007 


PG1115- 


h080 


2.0 


0.992 


0.991 


0.016 


0.015 


NOTE. - 


— Full table to 


appear 


in the on 


line edition. 



number of parameters involved in quadruple lenses, / qua d 
gives a rough criterion of how well our selection algorithm 
can identify quadruple lenses. 

5. REMOVING SINGLE QUASARS 

The nonzero completeness at 9 — > implies that our 
morphological selection algorithm falsely identifies "nor- 
mal" quasars as lenses. Indeed, when we apply our algo- 
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Fig. 14. — The distribution of inferred magnification diff erences (Am) and image separation (8) of simulated images, derived from 
fitting by two PSFs with GALFIT software IPeng eta l. 2002). The true image separation and the i-band flux ratio are set to 8 = 1" and 
fi = [10 — °' 5 , 1] (for double lenses), respectively. Crosses indicate double lenses, and open squares are quadruple lenses. The left and right 
panels show the results of u-band an d a- b and images, respectively. The dotted lines indicate our proposed limit in the image separation-flux 
ratio space (SI) defined in equation 1161 . 

sources by making a cut in image separation-flux ratio 
space. We do this fit on both the u- and z-band images: 
We use the z-band image because the flux ratio is defined 
by that of the z-band image. The u-band is used because 
it is a powerful discriminator between lensed pairs and 
superpositions of quasars and stars or quasars and galax- 
ies. 

To determine the region of parameter space we should 
exclude, we perform fitting of simulated images (see 
M4.ll) . We consider only the extreme image separation, 
9 = 1", for which fitting by two PSFs will be most diffi- 
cult. The flux ratio varied over the range fi = [10~ ' 5 , 1]. 
We show the results for both u- and z-band simulated im- 
ages in Figures 1141 (distribution in image separation- flux 
ratio space for each band) and !15l (distribution in u- and 
z-band image separation space) . It is seen that quadruple 
lenses show smaller image separations on average than 
double lenses; this is clearly because GALFIT tends to 
fit the two brightest images that are not necessarily the 
two separated farthest from each other. From this result, 
we determine the selection criterion to be the following: 

(SI): 9{u) > 0"4, 

6(i) > 074, 

Am(u) < 2.5, 

Am(i) < 2.2, 

\e(u)-e(i)\<i", (i6) 

where 9 and Am are the separation and magnitude dif- 
ference between fitted two PSF components. We select 
candidates that meet all of the above criteria. In our sim- 
ulation, more than 99% of double and quadruple lenses 
with 9 = 1" pass the condition, thus we conclude that 
this additional selection hardly changes the high com- 
pleteness of our statistical lens sample estimated in 2J 
The condition (SI) rejects ~ 85% of candidates selected 
by the morphological selection (see ^ ■ 

6. EFFICIENCY 







1 2 
0(u)[arcsec] 



and i-band image separations 
1" and the i-band flux ratio 



Fig. 15. — The distribution of 
(8) of simulated images with 8 

fi = [10 — °' 5 ,1] (for double l enses), derived fr om fitting by two 
PSFs with GALFIT software IPeng et al.M2002l) . Crosses indicate 
double lenses, and open squares are quadruple lenses. The dotted 
lines indicate our pr opos ed limit in the image separation space (SI) 
defined in equation 1161 . 

rithm to DR3 quasar catalog, we find many lens candi- 
dates, many of which appear to be point sources (see 
although most normal single quasars are removed by our 
morphological selection. Therefore, in this section, we 
present an additional selection which is aimed to exclude 
such point sources. 

To do so, we fit systems selected by our morphological 
selection algorithm with t wo PSFs using the GALFIT 
package l)Peng et al.l 12002}) . and derive the image sepa- 
ration and flux ratio. We use stars in the same field 
as a PSF template in fitting. The fitting of any single 
object with two PSFs should result in either too large 
magnitude differences or too small image separations. 
Here we only consider gravitational lenses included in 
the proposed statistical sample, i.e., those with 9 > 1" 
and fi > 10~ ' 5 . The point sources should have signifi- 
cantly smaller image separation and/or larger flux ratios 
than the above limit. Therefore we can reject such point 



A definitive estimate of the efficiency requires follow-up 
observations of all candidates which is still in progress. 
Nevertheless, in this section we discuss the efficiency of 
our selection algorithm in light of the follow-up observa- 
tions that we have performed thus far. 
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TABLE 4 

Preliminary efficiency of our selection algorithm 



6 N mor a 


iVmor + b 




N j d 

1 " cand 


N ohs c 


-^lens 


-^lens /-^cand 


< 2" 568 


88 


9 


96 


77 


8 


0.083 


2" - 7" 74 


8 


59 


66 


62 


2 


0.030 


7" - 20" 6 





155 


155 


141 


1 


0.006 


< 20" 648 


96 


223 


317 


280 


11 


0.035 



a Number of candidates selected by the morphological selection. 
b Number of candidates selected by the morphological selection plus additional 
selection presented in JF] 
c Number of candidates selected by the color selection. 

d Total number of candidates selected by the morphological and color selection. 
It does not necessarily coincide with the sum of 7V m0 r+ and N co \, because some 
of candidates are selected by both algorithms. 

c Number of candidates for which judgement of lens/nonlens was made on the 
basis of follow-up observations and/or various survey data. 

'Current number of confirmed lenses included among the candidates. 



To do so, we again adopt the DR3 quasar sample 

I) Schneider et all2005|) . As discussed in ® we use 22,682 
quasars at 0.6 < z < 2.2 and 15.0 < i COI < 19.1 as our 
source quasar sample. We apply our algorithms to the 
quasar sample, and select 648 candidates from the mor- 
phological selection and 228 candidates from the color se- 
lection. By applying the additional cut described in §^to 
morphologically selected candidates, the number of can- 
didates from the morphological selection reduces to 97, 
only ~ 15% of the original. For color selection, we only 
selected candidates with image separations smaller than 
20'.' 1 and magnitude differences smaller than 1.3, since we 
concentrate our attention to the lens sample in the range 
1" < 6 < 20" and f d > lO" ' 5 . It is expected that the ef- 
ficiency is a strong function of the image separation, thus 
we consider the three image separation ranges, 9 < 2", 
2" < 8 < 7", and 7" < 8, and derive the efficiency for 
each image separation range separately. The image sep- 
arations of morphologically selected candidates we used 
are those obtained by fitting i-band SDSS images with 
two PSFs using GALFIT (see 351- The result is summa- 
rized in TableQ] We note that the algorithm successfully 
identified all known lenses included in the subset of the 
DR3 quasar catalog, thereby confirming the high com- 
pleteness shown in 21 Since we still have candidates 
that are to be observed, the efficiency shown in Table 0] 
should in principle be regarded as the lower limit. How- 
ever, we believe that this is close to the true value because 
of the "observation bias" whereby we preferentially ob- 
serve better candidates. As expected, the efficiency is 
dependent on the image separation: The efficiencies are 
8%, 3%, and 0.6% for 8 < 2", 2" < 8 < 7", and 7" < 8, 
respectively. Most false positives, which have been re- 
vealed by follow-up obse rvations, were quasar- star pairs 
or quasar pairs (see also lHennawi et al"]l2006alK For the 
smallest image separation bin, the efficiency is compara- 
ble to that of the CLASS, which d iscovered 13 gravita - 
tional lenses out of 149 candidates l)Browne et alJ l20031. 
We note that the candidates can be restricted further by 
comparing with other existing observations (see Paper 

II) such as Faint Im ages of the Radio Sky at Twenty-cm 
l)Becker et al!l99 "5f) , which should increase the efficiency. 

7. CONCLUSION 

We have presented the SQLS lens candidate selection 
algorithm from the SDSS spectroscopic quasar catalog. 



Our algorithm consists of two selection methods, a mor- 
phological selection that identifies "extended" quasars, 
and a color selection that picks up adjacent objects with 
similar colors. Both selection methods rely only on imag- 
ing parameters that the SDSS reduction pipeline gener- 
ates, thereby allowing easy and fast selection of lens can- 
didates. To reduce the number of candidates further, we 
fit the images of morphologically selected candidates with 
two PSFs, deriving the image separations and magnitude 
differences, and applying cuts in the image separation- 
magnitude difference space. This allows us to remove 
most single quasars that are otherwise selected by the 
morphological selection. 

We have performed simulations of lensed quasar images 
in the SDSS to determine the effectiveness of the selec- 
tion algorithm. In the simulations, distributions of field 
parameters (i.e., seeings and sky levels) as well as those 
of source quasars are taken from the real SDSS data. 
We have also taken differential reddening by lens galax- 
ies into account. The simulated SDSS images are then 
passed through the imaging data processing pipeline used 
in the SDSS. This allows a reliable estimate of the com- 
pleteness of the algorithm to be made. We have found 
that our selection algorithm is almost complete at image 
separations larger than 1" and flux ratios larger than 
10~ - 5 . We have also quantified the magnification factor 
of lensed images as a function of the image separation, 
which is important for accurate computation of magni- 
fication bias. The algorithm successfully identifies both 
double and quadruple lenses, although there is a sys- 
tematic difference of completenesses between double and 
quadruple lenses. The efficiency depends strongly on the 
image separation; at separations less than 2" the prelim- 
inary efficiency of our selection algorithm is about 8%, 
comparable to that of the CLASS. 

An important effect we have neglected in our simula- 
tion is fluxes from lens galaxies. Although for quasars 
with i < 19.1 and redshifts above 0.6 lens galaxies 
tend to be fainter than lensed quasars, they could af- 
fect quantitative results. An important point is that 
lensed quasars are not selected by the quasar spectro- 
scopic target selection if the fluxes of lens galaxies dom- 
inate: A go od example is the hig h-redshift lens SDSS 
J0903+5028 (|Johnston et al.ll2003l) i that was targeted for 
SDSS spectroscopy as a luminous red galaxy, rather than 
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a quasar, because the lensing galaxy is brighter than the 
lensed quasar images. Therefore such systems will not 
be included in SDSS spectroscopic quasar catalogs. This 
means that we need to simulate quasar spectroscopic tar- 
get selection for lensed systems, as well as lens selection 
algorithms, to address the effect of lens galaxy fluxes. 
This is beyond the scope of this paper. 

The well-understood selection function, together with 
the knowledge of redshift and magnitude distributions of 
source quasars, is very important to control possible sys- 
tematic effects. We plan to construct a first SQLS com- 
plete lens sample from the DR3 quasar catalog, which 
will be presented in Paper II. 

We thank Issha Kayo for careful reading of the 
manuscript. N. I. is supported by JSPS through JSPS 
Research Fellowship for Young Scientists. M. A. S. is 
supported by NSF grant AST 03-07409. 

Funding for the SDSS and SDSS-II has been pro- 
vided by the Alfred P. Sloan Foundation, the Partic- 
ipating Institutions, the National Science Foundation, 
the U.S. Department of Energy, the National Aeronau- 
tics and Space Administration, the Japanese Monbuka- 
gakusho, the Max Planck Society, and the Higher Educa- 
tion Fundi ng Council fo r England. The SDSS Web Site 
is http:/ /www. sdss.org/ 

The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Par- 
ticipating Institutions are the American Museum of Nat- 
ural History, Astrophysical Institute Potsdam, Univer- 
sity of Basel, Cambridge University, Case Western Re- 
serve University, University of Chicago, Drexel Univer- 
sity, Fermilab, the Institute for Advanced Study, the 
Japan Participation Group, Johns Hopkins University, 
the Joint Institute for Nuclear Astrophysics, the Kavli 
Institute for Particle Astrophysics and Cosmology, the 
Korean Scientist Group, the Chinese Academy of Sci- 
ences (LAMOST), Los Alamos National Laboratory, the 
Max-Planck-Institute for Astronomy (MPA), the Max- 
Planck- Institute for Astrophysics (MPIA), New Mexico 
State University, Ohio State University, University of 
Pittsburgh, University of Portsmouth, Princeton Uni- 
versity, the United States Naval Observatory, and the 
University of Washington. 

APPENDIX 

OBJECT TYPES OF LENSED QUASAR 
SYSTEMS 

Gravitationally lensed quasars are sometimes classified 
as extended objects (objc_type = 3) by PHOTO. This 
is very important for lensing of high-redshift (z > 3) 
quasars because extended high-redshift quasars are not 
targeted by the quasar target selection pipeline (see Fig- 
ure [JJ. Therefore, we need to take this bias into account 
when we constrain the shape of the high-redshift quasar 
luminosity function from the absence of s trong gravita- 
tional lensing (e.g.. iRichards et aLl l2006a^. In addition, 
this bias strongly affects the search for lenses using a 
photo metrically selected quasar sample ( Ric hards et al J 
2004), because extended quasars are not included in the 
photometric sample, again in order to enhance the effi- 
ciency. In this appendix, we derive the fraction of lens 
systems that are classified as extended objects from our 



simulations of lensed quasar images (HJ. Although we 
focused our attention on bright (i < 19.1) low-rcdshift 
quasars and therefore the result may not be applied to 
high-redshift quasars that are fainter and have quite dif- 
ferent colors, the result is still useful to understand how 
PHOTO clas sifies lensed quasar images. 

Figure [ATI shows the fraction of double lenses that are 
classified as extended objects in the 9-fi plane. As ex- 
pected, the fraction is a strong function of the image sep- 
aration: The lenses with 9 ~ l'/2 are most likely judged 
as extended objects. The fraction depends on the flux 
ratio; lenses with larger /j are mo re li kely to be classi- 
fied as extended objects. In Figure lA"2l we plot the same 
fraction for quadruple lenses. Although the system- by- 
system difference is large, at 9 > l'/5 quadruple lenses 
are more likely to be classified as extended objects on 
average than double lenses, which is clearly because the 
existence of close fold/cusp image pairs. 
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Fig. Al. — Contours of the fraction of simulated double lenses 
that are classified as extended objects (objc_type = 3). The defi- 
nition of lines is same as in Figure Hfl 
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Fig. A2. — The fraction of quadruple lenses that are classified 
as extended objects (objc_type = 3). The dotted line shows the 
fraction of each quadruple lens image configuration and flux ratios 
taken from real observational data (see Table |21 . while the thick 
solid line indicates the fraction averaged over the 20 dotted lines. 
The fraction for double lenses, averaged over the flux ratios in the 
range 10~ 0,5 < /; < 1, is also shown by a dashed line for reference. 
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